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Abstract 

Ten-dimensional models, arising from a gravitational action which includes 
terms up to the fourth order in curvature tensor, are discussed. The spacetime 
consists of one time direction and two maximally symmetric subspaces, filled 
with matter in the form of an anisotropic fluid. Numerical integration of the 
cosmological field equations indicates that exponential, as well as power-law, 
solutions are possible. We carry out a dynamical study of the results in the 
Hext — Hint plane and confirm the existence of attractors in the evolution of the 
Universe. Those attracting points correspond to "extended" De Sitter space- 
times in which the external space exhibits inflationary expansion, while the 
internal one contracts. 

PACS Codes: 98.80.Hw , ll.lO.Kk 

I. Introduction 

The mathematical background for a non-linear Lagrangian theory of gravity was 
first formulated by Lovelock^, who proposed that the most general gravitational La- 
grangian is 

n/2 

^ = V-g ^{m)^{m) (1-1) 

where A(m) are coupling constants, n denotes the manifold's dimensions, g is the 
determinant of the metric tensor and are functions of the Riemann curvature 
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tensor, of the form 

r , , _ J_Xl3l-P2m 7?aia2 npOl2m-ia2m /-I qX 

Hm) — 2m"ai...a2m '^PiP^ '^/32m-i/32m y^-"^) 

where 6p is the Kronecker symbol, £(o) is the volume n-form which gives rise to the 
cosmological constant, £(i) = ^TZ is the Einstein-Hilbert (EH) Lagrangian and £(2) is 
the quadratic Gauss-Bonnett (GB) combination^. Euler variation of the gravitational 
action corresponding to Eq.(l.l) yields the most general symmetric and divergenceless 
tensor, which describes the propagation of the gravitational field and depends only 
on the metric and its first and second order derivatives^. 

While quadratic Lagrangians have been widely studied (e.g. see Refs. [3,4] and ref- 
erences therein), cubic and/or quartic Lagrangians only recently have been introduced 
in the discussion of cosmological models in the framework of superstring theories^"^*'. 
The reason is that, it is very hard to derive and (even harder) to solve the cor- 
responding field equations. In this case, solutions may be obtained only through 
certain numerical techniques^^'^^, where the idea of " attractor" plays a central role^^: 
If some special spacetime is the attractor for a wide range of initial conditions, such a 
spacetime is naturally realized asymptotically. Since the ten-dimensional superstring 
theory is a candidate for a realistic unified theory, it is very important to investigate 
whether a similar attractor exists in this theory. 

In the present paper we integrate numerically the field equations, resulting from a 
quartic gravitational Lagrangian, to obtain anisotropic, ten-dimensional cosmological 
models. The spacetime consists of one time direction and two maximally symmet- 
ric subspaces, FRW®FRW: The external space, representing the ordinary Universe 
and the internal one, constituted by the extra dimensions. The internal space is a 
compact manifold of very small "physical size" with respect to that of the "visible" 
space at the present epoch^^'^^. Since, on the other hand, at the origin the two sub- 
spaces were of comparable physical size, the internal one must have somehow been 
contracted towards a static value of the order of Planck length, Ipi ~ 10~^^cm, to 
achieve "spontaneous compactification"^^ . Compactification is a topological process 
of quantum origin, which leads to the separation of the extra dimensions from the or- 
dinary ones^^. In what follows we consider models of an already compactified internal 
space, i.e. we study only its contraction. 

In Section II we derive the explicit form of the field equations for a quartic theory 
in ten dimensions, in which both subspaces are filled with an anisotropic fiuid. In 
Section III we solve numerically the field equations, for a wide range of initial condi- 
tions and for several values of the "/ree" parameters involved, as regards (1) vacuum 
models of flat subspaces and (2) perfect fluid models of positively curved subspaces. 
Next, we carry out a dynamical study in the H^xt — Hint plane, where each Hj repre- 
sents the Hubble parameter of the corresponding subspace. Accordingly, we confirm 
the existence of attracting points and investigate their evolution with respect to the 
variation of the coupling constants X(m)- The explicit time-dependence of the un- 
known scale functions may be subsequently determined by solving the linearized field 
equations around these attracting points. The corresponding analysis is presented in 
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Section IV. 



II. The field equations in a quartic gravity theory 



We consider a ten-dimensional line element, representing cosmological models which 
consist of two homogeneous and isotropic factor spaces, of the form 



-dt^ + R^{t) 



1 + 



+ S\t) 



X/j=4 {dx'^ 



1 + \k 



int 



[2.1] 



where U — 1 — c, R{t) and S{t) are the cosmic scale functions of the external and 
the internal space respectively, k^xt = ~1,0, +1 is the curvature parameter of the 
"ordinary" space and kint = 0, +1 is the corresponding parameter of the internal 
one. Therefore, the extra dimensions may be compactified either in a six-dimensional 
sphere, for kint = +1, or in a six-dimensional torus, for kint — 0. The spatial section 
of the metric (2.1) can be viewed as the direct product of two FRW models with three 
and six dimensions respectively^. These models may be obtained through Hamilton's 
principle, from a ten-dimensional action in which the gravitational part is of the form 



I = 



1 



V, 



int 



'9 



A(o)>C(o) 



+ A(2)£(2) + ^(3)^(3) + ^{4)^(4) 



A (2.2) 



where each of C^m) is given by Eq.(1.2), A(^) are the corresponding coupling constants 
and Vint is a normalization constant^^, corresponding to the "physical size volume" of 
the internal space, once it may be considered static"^. The field equations read 



— 87rGin T, 



10 -Lnu 



(2.3) 



where is the Lovelock tensor up to the fourth order in curvature (Greek indices 
refer to the ten-dimensional spacetime)"*^ and Gio = G Vint is the ten-dimensional 
gravitational constant^^. T^,y is the energy- momentum tensor of an anisotropic per- 
fect fluid source, of the form T^^, = diag (p, —pext, —Pint, ■■■), where p is the total 
mass-energy density, while p^xt and Pint are the pressures associated to each factor 
space, separately. For the metric (2.1) Eq.(2.3) is decomposed into three independent 
equations of the form (cf. Ref. [9]) 



IGttGio P = A(o) + 6A(i) 
+ 72A(2) 



^ + 5g + 6(|)(|) 



+ 720A 



(3) 



+ 36g(|)^(|)^ 



+ 17280A(4) PQ 
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(2.4a) 



- 167rGio Pext = A(o) + 2A(i) 



P+15g + 12(|)(|) + 2(|) + 6(|) 



+ 24A 



(2) 



+ n^n^r + ^pq + 4og(|)(|) + 2og(|) 



+ 2P(-) + lOg(-) + 20(-)(-)(-) + 4(-^)(-^)(^) 



+ 720A(3) 
R 



g3 + 3Q^p + i2g^(|)(f ) + m^) + m^) 



+ i2g(^)^(|) + 4pg(|) + 8(f )(|)^(|)^ + 24g(|)(|)(f ) 



+ 8g(-)(-)(-) 



+ 5760A(4) 



2g^(f: 



+ i2g^(f )(f )(f ) + )(f )(f ) + 6^^'(f ) 
+ 24g(|)(|)^(|)^ + pg^ + 6g^(|)^(|)^ 



(2.46) 



167rGio Pint = A(o) + 2A(i) 



3P+10g + 15(|)(|) + 3(|) + 5(|) 



24A 



(2) 



+ 3og(|)(|) + iog(|) + 5P(|) + iog(|) + io(§)(§)(|) 



+ 720A 



(3) 



4(f)^(|)^ + 3pg^ + 3g^(|)(|) + 3g^(|) 



+ 2g^(|) + 6Pg(|)(|) + i2g(f )^(f )^ + 2pg(|) 

+ 6pg(|) + 4(f )(|)^(|)^ + i2(f )(f )^(f )^ + 4^(|)(f )(|) 



+ 5760A 



(4) 



6g^(f )(|)(f ) + 3pg^(|) + i2g(f )(|)^(f )^ 
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+ i2pg(|)(f )(|) + 3pg^(f )(|) + ^Qi^n^f 

+ 3Pg^(|) + 12g(|)(f + 8(§)(f (2.4c) 

where an overdot denotes derivative with respect to time and we have set 

^=(ff + ^ ' ^ = + ^ (2.5) 

Since the Lovelock tensor is divergenceless, C^^^^ = 0, we obtain the conservation law 
T^^ = 0, which gives 

R S 
p + 3{p + pext) + 6 (P + Pint) ^ = (2.6) 

Further inspection of the system of Eqs. (2.4) and (2.6) shows that only three of 
them are truly independent. Thus, the problem is completely determined by those, 
plus the two equations of state for the matter content, one for each subspace^''. In the 
present article we consider two cases with regard to the energy-momentum tensor: 

(a) Vacuum models, p = 0, in connection to flat spatial sections {kext = = kint) a-nd 

(b) Models of an hctcrotic supcrstring gas^", p^^t = \p and pint — 0, in connection to 
possitively curved spatial sections {kext — ^ — Knt) ■ In the later case, the conservation 
law (2.6) gives 

P = 1^ (2-7) 

where M is an integration constant. Thus, the external space is radiation dominated^'^^'^*^. 

In principle, wc may integrate the system of Eqs. (2.4) and (2.6) to obtain the 
form of the unknown scale functions. However this is not an easy task, even in the 
most simple and symmetric cases^. Nevertheless, we may get a good estimation of 
their dynamic behaviour through numerical integration^^. 

Once the two equations of state are determined, Eq.(2.6) may be readily solved to 
give the unknown energy density and pressures, as functions of R{t) and S{t). These 
expressions are subsequently introduced in the r.h.s. of Eqs. (2.4). Now, only two of 
these equations are truly independent. The third one corresponds to a constraint, 
to be satisfied by the solutions of the system. As such, we choose Eq.(2.4a). The 
remaining independent field equations (2.46) and (2.4c) may be recast in the form of 
a first order system (see also [11]), as follows 

Hext = Gi {Hext 1 Hint ■, X , Y) (2.8a) 

Hint = G2 {Hext , Hint , ^ , Y) (2.86) 

X^-XHe,t (2.8c) 
Y^-Y Hint (2.8d) 
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where we have set 



next — , -"mt ~ ~g ' ' ~ ^52" ^ 

and the exphcit forms of the functions Gi and G2 are given in the Appendix A. 

Finally, it is convenient to make a parameter rescaling in the field equations, of 
the form 

«^ = ^ , m = 0, 1, 2, 3, 4 (2.10) 

A(i) 

where A(i) = (IGttG)"^ is the coupling constant in the four-dimensional General 
Relativity (GR). The value of the normalized coupling constants. Urn {i^m ^ 1); is 
directly proportional to the contribution in the field equations of the corresponding 
m~th order non-linear term, with respect to the results obtained in the EH cosmology. 
Clearly, ki — 1. 



III. Numerical Results 

We integrate numerically the system of Eqs.(2.8). The constraint (2.4a) is checked 
to be satisfied with an accuracy of 10~^° along integration. The initial conditions 
i^Q^*, i?o"S -'^o, ^0 are chosen so that: (a) Xq — Yq, i.e. at the origin, the two factor 
spaces are separated, but of the same "physical sizc"^^. (b) H^^^ > 0, i.e. initially the 
ordinary space expands, in accordance to what we observe at the present epoch^^'^^. 
(c) Hq"'^ < 0, i.e. at the origin, the internal space contracts, in correspondance to 
" spontaneous compactification" 12,16,17,21 The cases where either i^Q^* < Oori^o"* > 
are not permitted, since the constraint equation is not satisfied. Nevertheless, the case 
where both conditions H^^^ < and ifg"* > are valid is acceptable by numerical 
analysis. Actually, it corresponds to the time-reversed solution of the system (2.8). 

The time coordinate is measured in dimensionless units, being normalized with 
respect to the Planck time, t — t/tpi {tpi = \fG ~ 10~^^sec). The limits of numerical 
integration range from r = to r = 10^. The upper limit coincides with the origin 
of the GUT epoch^^, Iqut = 10^ tpi, corresponding to the end of the string regime^^. 
However we have to point out that, although the origin of the time coordinate is set 
at r = 0, the equations (2.8) may not be valid in the region < r < 1 since, in the 
absence of a quantum gravity theory, there is always a region of ambiguity around 
t — 0, of the order of Planck time^^~^^. 

The solution of the system (2.8) may be represented as curves in the /fp,,.^ — Hint 
plane. Any point located on these curves always satisfies the constraint condition 
(2.4a). Thus, the curves actually represent "orbits" of the dynamical system under 
study. Each curve, corresponding to a different set of initial conditions, is bounded by 
fixed points (or infinities) and represents a different type of evolution for the Universe. 

In what follows, we focus attention on the existence and the evolution of attracting 
points in the Hg^t — Hint plane. The reason rests in the physical meaning of the 
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attractor: No matter what the behaviour of a cosmological model at the origin might 
be, it will always end up to evolve as indicated by the location of the attracting point 
in the H^xt - Hi^t plane. 

(1) Vacuum models with spatially flat subspaces 

We study the evolution of vacuum ten-dimensional cosmological models, with metric 
of the form (2.1), in which both subspaces are spatially flat, i.e. k^xt = = hint- 

Thus, X = = y. 

The first case to study are the GB models (see also Refs. [11,12]). In this case, 
K2 — 1 and Ko = = /t3 = /t4. The non-linear curvature contributions to the field 
equations come out from the quadratic terms alone. The time evolution of the Hubble 
parameters is presented in Fig. la. We see that both parameters evolve to approach 
constant values in the later stages. This situation verifies the existence of attracting 
points in the Hg^t — H^nt plane during the evolution of the Universe. Therefore, for 
a wide range of initial conditions, both subspaces will end up to evolve as De Sitter 
spaces, in complete correspondence to the results of Ishihara^^. 

We also observe that Hext > and Hint < 0. Therefore, while the internal space 
contracts exponentially to achieve spontaneous compactification, the external one 
expands, a fact that corresponds to an infiationary phase. This result indicates that 
the introduction of the non-linear curvature terms into the gravitational action may 
play an important role as far as the inflation is conserncd^^^'^°. The explicit location 
of the attracting point is shown in Fig. lb. The attractor corresponds to the flxed 
point D2 recognized by Ishihara^^ in the evolution of the extended De Sitter models 
in GB theory. 

The next step is to introduce into the problem a "bare" cosmological constant, 
A, corresponding to the expectation value of the vacuum energy density^^. Now, in 
addition to /t2, we also have 7^ 0, while k,^ = = K4. When Kq E [0 , 1] the value 
of the cosmological constant in physical units is A = 2k,o x 10^"^^ cm~^, which is quite 
small. 

The behaviour of the model is qualitatively similar to the previous case. Again we 
verify the existence of an "attractor". Both subspaces correspond to De Sitter models. 
The external space exhibits inflationary expansion, while the internal one contracts. 
However, in this case, the location of the attracting point D2 has changed to higher 
absolute values in the evolution of H^xt and Hint (Fig- 2a). We may determine 
explicitly the law of the attractor's displacement in the Hext — Hint plane, caused by 
variations of the cosmological constant. 

In general, to determine the exact location of the attracting points in an Hf.xt—Hint 
plane, requires to set 

G^{H,xt,Hint,X,Y)^^ (3.1.1a) 
G2{Hext,H,nt,X,Y) = Q (3.1.16) 
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In the case of flat and vacuum subspaces [X — Y — p^xt — Pint — 0), Eqs.(3.1.1) read 



fl{Hext^Hint^K,rn) — [Gl2G2Q — G22GlQ\x=Y=Q ~ ^ (3.1.2a) 
f2{HexuHint,K,m) — [G21G1Q — GiiG2o\x=Y=0 — (3.1.26) 

where m = 0,2,3,4 and the quantities Gij are presented in the Appendix A. We 
differentiate the functions /i and /2 with respect to Hext , Hint and Km, to obtain a 
system of first order differential equations {"variational equations") 

dfi 9fi 9fi 

= i arr^ )pdHext + ( orr^ )pdHint + y^{^^)pdKm = (3.1.3a) 



df2 , , r df2 , , ^,df, 



df2 = {i^)pdH,,t + {^)pdHi^t + Y,{^)pdKm = (3.1.36) 



The system (3.1.3) may be used, to determine the evolution of the attracting point 
-D2 {Hext, Hint), under the variation of the normalized coupling constants Km- For 
K2 — 1, in the case of vanishing k^ and K4, the evolution of the attractor D2 {Hext, Hint) 
with respect to the variation of the cosmological constant kq, is given by 

(^) - (3.1.4a) 

(^) . (^) (3.1.46) 

where the functions PP, QQi and QQ2 are given in the Appendix B. Subsequently, 
the system (3.1.4) is evaluated by numerical integration. The corresponding results 
are shown in Fig. 2b. Using least square fitting, we see that the displacement of D2 
takes place along the straight line 

Hint = -0.075Hext - 0.071 (3.1.5) 

The investigation of the behaviour of the models under consideration by including 
a third order curvature term, corresponds to study them at earlier epochs in the 
history of the Universe. Indeed, if we are interested in the behaviour of the model 
very close to the initial singularity, the leading terms to consider in the field equations 
are those with the highest power in Q^, i.e. those obtained from the highest order 
terms in the gravitational action^. 

The time-evolution of the model is quite similar to the previous cases. In the later 
stages it corresponds to an extended De Sitter model, in which both subspaces evolve 
exponentially. The external space expands, while the internal one contracts (Fig. 3a). 

Again, we verify the existence of an attracting point P in the evolution of the 
Hubble parameters and we investigate its behaviour as k^ increases, from to 1, i.e. 
until it becomes as important as the quadratic term. The evolution of the attractor in 
the >C(3)-theory, with respect to the variation of K3, may be obtained in a similar way 
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as in the kq case. We differentiate tlie functions /i and /2 witli respect to Hext , Hint 
and to obtain a first order system of differential equations wliicli, for K2 = 1 and 
for vanisliing Kq and K4, will determine the displacement of P in the Hgr^f — Hint plane, 
under the variation of ^3. 

The corresponding results are presented in Fig. 3b. We observe that the attractor 
moves to higher absolute values of Hint as increases. This result has a clear physical 
meaning. Since increasing K3 corresponds to study the earlier stages in the evolution 
of the Universe, we see that at these epochs the internal space contracts at higher 
rates than those of the GB theory. Then Fig. 3b verifies that at the late stages, where 
the GB theory holds alone, the value of the internal Hubble parameter decreases in 
order to achieve stabilization. 

Again, the law of displacement of P in the H^xt — Hint plane may be estimated 
using best- fit methods. In this context, we find that it may be represented by a sixth- 
order polynomial H^xt = PeiHint), with coefficients: Gq = 0.7373 , Oi = —25.594 , 02 = 
-457.324 , aa = -3323.9 , 04 = -12156.9 , 05 = -22135.2 and oe = -15905.6. 

Finally, to solve the cosmological field equations when all terms in the action (2.2) 
are included (i.e. K4 7^ 0), corresponds to study the dynamic behaviour of the model 
under consideration at even earlier epochs. The results are slightly different from 
those of the previous case (Fig. 4a). Again, in the later stages, the model consists 
of two De Sitter subspaces and there exists an attracting point. The attractor 's 
displacement in the H^xt — Hint plane is obtained in a way similar to the and 
cases and may be represented by a third-order polynomial. Hint = PsiHext), with 
coefficients: = 7.47, 61 = —34.27, 62 = 51.33 and 63 = —25.74. The corresponding 
result is shown in Fig. 4b. 

Hence, we may conclude that in every case where non linear terms are included, 
the "extended" De Sitter solution (i.e. an exponentially expanding external space in 
connection to an exponentially contracting internal one) corresponds to an "attractor" 
of the dynamical system under consideration. Accordingly, (in our model) no matter 
how the Universe may originate, there is at least one period during its time-evolution 
in which it exhibits inflation of the ordinary space, accompanied by spontaneous 
compactiflcation of the internal one^^'^°. 

(2) Perfect fluid models of curved subspaces 

We consider a ten-dimensional metric of the form (2.1), which now represents a class 
of cosmological models with positively curved subspaces {kext = 1 = kint)- Then, 
X = R^^{t) and Y = S^^{t) and we study the time-evolution of the cosmological 
models as results from the solution of the system (2.8). 

The numerical analysis is carried out in the same fashion as in the previous case 
of vacuum models. We consider that at the origin both subspaces are of the same 
"physical size", i.e. Xq — Yq, but they have different expansion rates, i?Q^* and ifg"*. 
As such we choose the corresponding range used in the vacuum case. We normalize 
both scale functions R{t) and S{t) to unity, with respect to their value at the Planck 
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epoch. That is 

R(t) ^ M , S(t) ^ M (3.2.1) 

where i?p; = Spi. As initial conditions we choose i?o = 100 = 5*0. 

We represent the matter fiUing the Universe by a closed or heterotic superstring 
perfect gas, with the following equation of state, deduced by Matsuo^° 

Pext ^ \ p ^ Pint = (3.2.2) 

Thus, the external space is radiation-dominated, while the internal one is pressure- 
less. It has been recently shown that, in this case, the two subspaces are completely 
disjoint^'^^. The time-evolution of the total mass-energy density p is accordingly given 
by Eq.(2.7). 

As regards the GB models {k,^ = = k,4^), we have performed a number of com- 
putational runs, varying the initial values of the Hubble parameters and the coupling 
constant k,2 as well, from /«2 = 0.1 to /«2 = 1- Numerical results in this case indicate 
that there is a considerable difference with respect to the vacuum-fiat models. It rests 
in the fact that the range of values of the coupling constant k,2 may be splitted into 
two parts. Each one of these parts leads to a different time-evolution of both the 
external and the internal scale functions. 

The first part consists of values of k,2 in the interval 0.1 < K2 < 0.65, i.e. when 
the contribution of the quadratic curvature terms is relatively small. In this case we 
expect that the time-evolution of the Universe will be only slightly different from the 
corresponding EH one. Indeed, the numerical results indicate that the system (2.8) 
admits solutions with a power law dependence of the scale functions upon time, of 
the form 

R{t) oc t""' (3.2.3a) 

S{t) oc r""^ (3.2.36) 

where the values of the indices mi and 1712 arc continuously increasing in the ranges 
0.25 < mi < 0.55 and 0.01 < m2 < 0.11, as K2 increases from 0.1 to 0.65. In 
this case, there are no attracting points in the evolution of the Universe. The last 
values in those ranges (0.55 and 0.11, respectively), both corresponding to the value 
K2 = 0.65, represent a Kasner-type regime^^'^^"^^ of the GB models. Indeed, the 
analytic approach in this case suggests that the two subspaces evolve as 

R{t) ~ tP' , S{t) ~ t-P^ (3.2.4) 

where both pi and p2 are possitive and in a ten-dimesional spacetime they satisfy the 
conditions 

3pi - 6p2 = 1 (3.2.5a) 
3pl + 6p^ = 1 (3.2.56) 
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The only physically acceptable solution of the system (3.2.5), compatible with the 
condition pi , P2 > 0, is 



Therefore, when K2 — 0.65, although the spatial sections are curved, the time- 
evolution of the Universe admits a Kasncr-typc solution. This solution actually lies 
on the interface between two different types of cosmological behaviour (Figs. 5a and 



The second type of time-evolution arises when 0.65 < K2 <!■ Then the Universe 
behaves, again, as an extended De Sitter spacetime (where the external space expands 
while the internal one contracts, both exponentially). In this case there exists an 
attracting point as in the vacuum-flat models (Fig. 6a). 

In conclusion, for a curved ten-dimensional GB cosmological model, filled with 
matter in the form of a superstring perfect gas, we may obtain three different types 
of cosmological behaviour, depending on the exact value of the normalized coupling 
constant K2: 

(a) Power-law solutions, with no attracting points, when 0.1 < ^2 < 0.65. 

(b) A Kasner-type model, when K2 — 0.65. 

(c) Extended De Sitter models, with an attracting point, when 0.65 < < 1- 

In all cases, the external space expands, while the internal one contracts. The 
inclusion of the contribution of the third and/or the fourth order terms in the field 
equations simply amounts to a modulation of those results (Fig. 6b). 



Analytic expressions, for the time-evolution of the model Universe considered, may 
be obtained by solving the cosmological field equations (2.8) around the attracting 
points. Accordingly, we investigate the cosmological behaviour of a vacuum, ten- 
dimensional model with spatially flat subspaces {X — Q — Y) within the context 
of the quartic Lagrangian theory under consideration. Clearly, setting some of the 
coupling constants A(„j) equal to zero corresponds to reducing the general theory to 
its lower case counterparts (EH- cosmology, GB-theory etc.). 

Since we are interested in the behaviour of the model around the attracting points, 
we consider the linearized equations 



pi = \ = 0.555 , P2 = ^ = 0.111 



5b). 



IV. Analytic Results 



ext 



(4.1a) 



Hint = ^ + H2{t) 



(4.16) 
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where Ai and A2 are the coordinates of the attractor, while Hi(t) and H2{t) represent 
small perturbations around those values {\Hi\ , \H2\ <^ 1). Therefore, to obtain the 
time-evolution of H^xt and H^nt, we only have to solve the system (2.8) linearized with 
respect to Hi and H2. 

The system of the cosmological field equations (2.8), hnearized with respect to Hi 
and H2, may be written in the form 

(3iHi + (32H2 + (33 



Hi{t) = 



aiHi + a2i?2 + ^3 

llHi + 72-^2 + 73 



(4.2a) 
(4.26) 



aiHi + Q;2-f^2 + «3 

where ctj,(3j and 7^ {j = 1,2,3) are constants, calculated directly from the lineariza- 
tion of the original equations, which depend on Ai , A2 and X(^m) (m = 0, 1,2,3,4). 
Prom Eqs.(4.2) we obtain 

dHi ^ /3iHi + 32H2 + ds 

dH2 71-H'i 72-ff2 + 73 
The solution of Eq.(4.3), in connection to Eqs.(4.1), will give, in the linear approxi- 
mation, the analytic expression of H^xt in terms of Hint- To solve Eq.(4.3), we need 
to have the solution {hi , h2) of the algebraic system 

piHi + P2H2 + p3 = (4.4a) 

jiHi + J2H2 + 73 = (4.46) 



We choose 



and furthermore, we set 



(4.5) 



w 



Hi - hi (4.6a) 

z = H2-h2 (4.66) 

We verify that the solution of Eq.(4.3) depends on several algebraic combinations of 
the constants aj , and jj, something that leads to several conditions between the 
coupling constants X(m)- Therefore, we consider the following cases: 

(a) /3i -|- 72 7^ 0: This combination corresponds to the most general case. Setting 



the solution of Eq.(4.3) reads^^ 



4/?27i + iPi - 72)' 



f /3i+72 l„ ,gl-72-^/^-27l^ 
/3l-72+^/^-27l^ 



(4.7) 



In 



/32^^ + iPl - ^2)ZW - -fiW" 



2(/3i+72) 



(/3l-72)-27l- 



for A < 
for A = 



2toarctan^^i^^^^ for A > 



(4.8) 
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where c is an arbitrary integration constant. 

(b) /3i + 72 = and A = with /327i < 0: In this case wc may proceed to derive 
the exphcit time-dependence of the Hubble parameters and the corresponding scale 
functions for both subspaces. Prom Eq.(4.8) we obtain 

H2 = CiHi + C2 (4.9) 

where 

Cl = ^^, C2 = 1 + /l2 - Ci/li (4.10) 

Now, Eq.(4.9) is inserted into Eq.(4.2a) to give 

where the constants S , e , ( and rj stand for the combinations 

S = Pici + (32 , (33 + (3iC2 

C = aiCi + a2 , T] = as + a\C2 (4-12) 
We consider the following cases: 

(i) 5 , C 7^ 0: In this case, Eq.(4.11) results in 

where to is an integration constant. Now, Eq.(4.1a) in connection with Eq.(4.13), 
may be easily integrated to give the form of i?(t), when the condition 

^ - ^ = (4.14) 

holds. Then, we obtain 

\nR{t)r. A,{t-to) + —{t-tof (4.15) 
which introduces a quadratic correction to the expected De Sitter solution. 

(ii) 5 , 77 7^ and C = 0: In this case we rediscover the solutions of Ishihara^^, 
obtained in the GB theory, as a particular case of the general solution. Indeed, from 
Eq.(4.11) we obtain 

Hi ^ Ce^ - ^ (4.16) 
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where C is an arbitrary integration constant. Therefore the corresponding external 
scale function is of the form 



In R{t) ~ - - to) + C|ef (4.17) 
Since the external space expands, we must have Ai > |. For C — 0, Eq.(4.17) reads 

R{t) - e(^^-t)(*-*°) (4.18a) 
corresponding again to a De Sitter phase, while for ^ <S 1 it yields 

R{t) ~ e(^i-t)(*-*°) (l + C(|)2ef (4.186) 

For e = Eq.(4.18b) corresponds to the solution of Ishihara (Eq.(15) of Ref. [12]) 
obtained in the framework of the GB theory. 

(iii) S — , C 7^ 0- Finally, in this case, Eqs. (4.1a) and (4.11) result in 



V 



3/2 



In R{t) ~ l^^i -fj{t- to) ±^[v' + 2eC(i - to)\ (4.19) 

where, in connection with the numerical results we must have Ai > ^. 

In concluding, we see that the coupling constants X^rn) may not be arbitrary. In 
every case, they should satisfy certain algebraic relations, depending on the form of 
the corresponding solution around the attracting points. 

Since both Eqs. (4.2) are almost of the same functional form, in all of the preceding 
cases, similar functional results may be obtained for the internal space, through the 
solution of Eq.(4.2b). In this case, however, we must take into account the fact that 
the numerical results indicate that the extra dimensions contract {A2 < 0). This 
argument may lead to additional constraints on the coupling constants A(^). 



V. Discussion and Conclusions 

In the present paper we have studied the time evolution of anisotropic, ten-dimensional 
cosmological models in the framework of a quartic Lovelock-Lagrangian theory of 
gravity^'^"^^. The cosmological models under consideration consist of one time di- 
rection and two homogeneous and isotropic subspaces: A three-dimensional external 
space, which represents the ordinary Universe, and a compact internal space, which is 
constituted by the extra dimensions. The evolution of the Universe depends on four 
free parameters. These are the coefficients \rn) which introduce the extra curvature 
terms in the gravitational Lagrangian (m = 0, 2, 3, 4). They are to be regarded as the 
coupling constants^'^. Since we have considered models of an already compactified 
internal space^^, we accordingly examine the process of its contraction^^"^^'^^. 
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The Universe is filled with matter in the form of an anisotropic perfect fluid. Given 
an equation of state for the matter content of each subspace, the time evolution of 
the Universe is completely determined by a system of three second-order, non-linear 
differential equations, consisting of the field equations (2.46) and (2.4c) together with 
the conservation law (2.6). The initial value field equation (2.4a) corresponds to a 
constraint which should be satisfied by the cosmological solutions. As regards the 
energy momentum tensor, we have considered two cases: (a) Vacuum models, p = 0, 
in connection with spatially flat subspaces and (b) Models of an heterotic superstring 
gas, with pext = |p and pint = 0, in connection with positively curved subspaces. 

The three independent equations may be subsequently expressed in the form of 
a first order system, Eqs.(2.8), involving the Hubble parameters Hext , Hint and the 
corresponding scale functions R(t) , S{t) of the two factor spaces. This system is 
evaluated numerically, for a wide range of initial conditions of the form Hq^^ > and 
jjmt ^ g solutions may be represented by curves in the H^xt — H^nt plane. Those 
curves correspond to the "orbits" of the dynamical system under study and each one 
of them, associated with a different set of initial conditions, represents a different 
type of evolution for the Universe. 

In the case of vacuum models with flat subspaces {k^xt = = /cj„t), the numerical 
results indicate that for all values of the coupling constants involved and also for a 
wide range of initial conditions, the Universe will always end up to evolve according 
to an extended De Sitter solution, i.e. an exponentially expanding external space, 
accompanied by an exponentially contracting internal one. Indeed, in this case the 
Hubble parameters of both subspaces approach constant values in the later stages. 
We have conflrmed that those values actually represent the attracting points of the 
dynamical system under consideration^^"^^. The appearence of attractors in the so- 
lution of the cosmological field equations is very important, since, if a spacetime is an 
attractor for a wide range of initial conditions, then it may be realized asymptotically 
in the later stages^^'^"^. Those results indicate that the existence of the non- linear 
curvature terms in the gravitational action may lead to inflation without the use of 
any phase transition^^'^^"^^'^^. 

Furthermore, we have investigated the evolution of the attractors under the varia- 
tion of the normalized coupling constants Km = ^{m)/^{i) (m = , 3 , 4). In all cases, 
the attracting points arc displaced at higher absolute values of Hint as increases 
from to 1. As regards the variation of k.^ and k.^, this result has a clear physical 
meaning. 

In determining the cosmological behaviour of the model very close to the initial 
singularity, the leading terms to consider are those with the highest power in (j), i.e. 
those obtained from the highest order terms in the gravitational action^'^. Therefore, 
the increase of (m = 3 , 4) corresponds to a more accurate study of the earlier 
stages in the evolution of the Universe^. Then, from Figs. 3b and 4b, we see that at 
those epochs the internal space contracts at higher rates than those of the GB theory. 
This means that in the later stages of the time evolution, where the GB theory holds 
alone, the absolute value of the internal Hubble parameter decreases, in order for the 
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extra space to achieve stabilization at a small physical size^^~^^'^^ '^^. 

In the vacuum case it is possible to derive the analytic dependence of the scale 
functions upon time, by linearizing the field equations around the values of i^ext and 
Hint at the attracting points. The corresponding results indicate that the functional 
form of the analytic solution depends on several algebraic conditions between the 
coupling constants A(ot) . Therefore, in a Lovelock-Lagrangian theory of gravity, the 
coupling constants may play an important role in determining the cosmological be- 
haviour of the model Universe. Nevertheless, the coefficients of each term in the 
Lagrangian either should be determined experimentaly or they should be given by 
some underlying foundamental theory^. In this context, we have rediscovered the 
solutions of Ishihara^^, obtained in the framework of the GB theory, as particular 
solutions of the general quartic case. 

The cosmological models with matter in the form of a superstring perfect gas, in 
which both subspaces are possitively curved {kext = 1 = Knt)-! can be treated only 
numerically. In this case, the evolution of the GB models depends additionally on the 
exact value of the normalized coupling constant K2. We have obtained three different 
types of cosmological behaviour: 

(a) Power-law solutions, with no attracting points, when 0.1 < fi;2 < 0.65. 

(b) A Kasner-type model, when K2 — 0.65. 

(c) Extended De Sitter models, with an attracting point, when 0.65 < ^2 < 1- 

In all cases, the external space expands, while the internal one contracts. The 
inclusion of the contributions of the third and/or the fourth order terms in the field 
equations, simply amounts in a modulation of the above results (e.g. see Fig. 6b). 
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The cosmological field equations (2.46) and (2.4c) may be recast inthe form of a first 
order system, as follows 



Appendix A 



ext 



Gl {Hext 1 Hint 1 X 1 Y) 



{A.la) 




{A.lh) 
{A.lc) 



X^-XH, 



ext 
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Y^-Y Hint (A.ld) 

where 

_ ^ext -^2 _ ^int t A 0\ 

The functions Gi and G2 are given by the expressions 

^ _ [Gi2{16lTPint + G20) - G22{16'n-pext + Gio)] ^ A o \ 

— fTw^ ^ ^ 1 (A.daj 

^ [G2l{l6'KPext + <^io) - GiillGTrp^^t + G'20)] 
[(J11G22 ~ G12G21] 

where we have set 

do = Ao + 2Ai5io + 24A2520 + 72OA3530 + 576OA4S40 (AAa) 

Gn = 4Ai + 24A2521 + 72OA3531 + 576OA4S41 (AAb) 

Gu = 12Ai + 24A2S22 + 72OA3S32 + 576OA4S42 (AAc) 

G20 = Ao + 2A1C10 + 24A2C20 + 72OA3C30 + 576OA4C40 (AAd) 

G21 = 6A1 + 24A2C21 + 72OA3C31 + 576OA4C41 (AAe) 

G22 = lOAi + 24A2C22 + 72OA3C32 + 576OA4C42 (AAf) 
and the quantities Bij and Cy are given by 

Sio = 3Hlt + 21Hl, + 12HextHint + X'' + 15Y^ (A5a) 

B20 = 35iff„i + 27Hl^^Hf^^ + mHextHfnt + '^Hl^tHint 
+ 7i/f^tX2 + 5y2 j^3y2 ^ X^) 

+ Y\bml, + lbHl^ + Am,,,H,nt) {AM) 

Bso = 7iff„t + ?,?,Hl^^Hf^^ + SH^^^Hf^^ + SeH^xtHf^t 
+ 9Y' {Hi, + Hi,) + SX'Y' + y« 

+ (l5-fff„i + 3AHl,Hl, + ASHextHlt + ^HltHint) 

+ 10i/2^,X2r2 + 7i/4,X2 + UHextHintY' (A.5c) 

S40 = 39HI,HI, + l2HextHlt + l2Hl,Hl, 



2 I \-2t^ 

3 I 1 o 0-3 



+ 3Hl,Y^ + X'Y 



+ Y [21H^^,H^^, + 12HextHint + '^'^H^^^Hi^t 
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int-' 



B21 = 10Hl^ + 4He,tHint+10Y^ 



+ Y^(l2Hl, + 8HeMnt) 



B41 = 2H-^^ + 12i/ext-f^i„t 

+ Y' (6//4, + 2AH,,tHl,) + 2Y' 



B22 — 20i?j^„j + 2ifg^j + 20HextHint 
+ 2X^ + 20Y^ 



B32 = GHfnt + ^"^H^xtHfut + '^^HextHf^t 

+ Y^(l2H,xtH,nt + QHl?j 
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4 
ex 



a 



40 



+ 3y^ (SHl^Hint + H^^H^^^ + H^) 
+ SX^y^ (iff„4 + Hl^ + HextHint) 

+ + 5HextH-^t + HltHtnt) 
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int 



C31 = 3i?4t + 6Hl^Hf^^ + 12HextHi^t 
+ 2Y' [sHl, + Hi, + 6H,xtHint) 

+ 2X^Hf^, + {3Y^ + 2X^Y^) 



a 



41 



+ 3F^ (2HextHint + Hi) 

+ 6Y'{2HextHl, + 3Hl,Hl,) 
+ QHl.X^Y^ + ZHLX^ + 3XV 



32 



= 2if4t + ISHl^Hl, + AHl^Hint + UHextHf^t 
+ 2Y' {2Hl, + 3HI, + 6i/e.t/^int) 

+ 2X^ (3Hl, + 2HextHint) + 2Y^ + 6X^Y^ 



19 



int 



Appendix B 



where we have set 



For fi;2 = 1, in the case of vanishing and ^4, the evolution of the attracting point 
D2 with respect to the variation of the cosmological constant kq, is given by 

The functions PP, QQi, and (5(52 are given by the expressions 

PP = P11P22 - P12P21 (P.2a) 

QQi = P1P22 - i^2i^i2 {B.2b) 
QQ2 = P11P2 - P1P21 (P.2c) 

Pi = -[2 + 24(10i/i^, - 3HlJ] (B.Sa) 
P2 = -[2 + 24(i/2^, + GHeMnt)] {BM) 

(P.3c) 
(P.3d) 
(P.Se) 
{B.3f) 

Now, the quantities Gij are given by 

G'n = 4 + 24(P2i)x=y=o (BAa) 

G12 = 12 + 24(P22)x=y=o (5.46) 

do = + 2(Pio)x=y=o + 24(P2o)x=y=o (5.4c) 

G'2i = 6 + 24(C72i)x=y=o (5-4o?) 

^22 = 10 + 24((722)x=y=o (5-4e) 



i^ll 


— G121G2Q 


+ ^12^201 


— (^221^10 


— G22G1QI 


P12 


— ^122^20 


+ G12G2Q2 


— (^222^10 


— (j 22^102 


P21 


= G2I1G1Q 


+ (^2lGioi 


— (^111(^20 


— (^11(^201 


-^22 


= Cj212G'lO 


+ ^21(^^102 


— (j^ll2(j^20 


— (j^ll(j^202 
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G20 = «o + 2(Cio)x=y=o + 24(C2o)x=y=o (BAf) 

where {Bij)x=Y=o and (Cjj)x=y=o denote the form of the corresponding quantities 
for X = = y and the symbols Gijk stand for 

Gijk = {-gj^)x=Y=o (-B.5) 
in which, /c = 1 for Hg^t and A; = 2 for i^^nt- 
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Figure Captions 



Fig. la: The evolution of tlie Hubble parameters in the GB theory [k2 = 1), for three 
different sets of initial conditions {H^^^^W^^. (0.75,-0.25) [solid line] (0.85,-0.15) 
[dashed hne] and (0.95, —0.05) [squares]. The time coordinate is measured in units of 

10^ tpi. 

Fig. lb: The orbits (in the Hg^t—Hint plane) of the dynamical system determined 
by the cosmological field equations for a model with flat subspaces, for three different 
sets of initial conditions. All orbits end at the attracting point D2 (0.8866, —0.1375). 

Fig. 2a: The change in the location of the attractor D2 when a non-zero cosmo- 
logical constant is included, for kq ~ 0.5. 

Fig. 2b: The displacement of the attractor on the Hext — Hint plane for a 
wide range of values of the cosmological constant (squares). Notice the very good 
agreement with the least square fitting result Hint = —0-075Hext — 0.071. 

Fig. 3a: The evolution of the Hubble parameters in the £(3)-theory, for K3 — 
0.15 and for three different sets of initial conditions (ifg^*, if™*): (0.75, —0.25) [solid 
line] (0.85,-0.15) [dashed line] and (0.95,-0.05) [squares]. The time coordinate is 
measured in units of 10* tpi. 

Fig. 3b: The displacement of the attractor on the H^xt — Hint plane for a wide 
range of values of the coupling constant K3 (squares) . There is a very good agreement 
with the plot of a sixth-order polynomial of the form H^xt = {Hint)- 

Fig. 4a: The evolution of the external Hubble parameter for Hq^^ = 0.85, when 
several combinations of the non-linear terms are gradually included in the field equa- 
tions. The time coordinate is measured in units of lOHpi. 

Fig. 4b: The displacement of the attractor on the H^xt — Hint plane for a wide 
range of values of the coupling constant K4 (squares). Notice the very good agreement 
with the plot of a third-order polynomial (solid line) of the form Hint = PsiHext)- 

Fig. 5a: The time-evolution of the positively curved external space, for several 
values of the normalized coupling constant K2- 

Fig. 5b: The corresponding evolution of the positively curved internal space. In 
both figures the time coordinate is measured in units of 10^ tpi. Notice the change in 
the cosmological behaviour of both subspaces when K2 < 0.65 and when K2 > 0.65. 

Fig. 6a: The orbits (in the H^xt — Hint plane) of the dynamical system determined 
by the cosmological field equations for a model with positively curved subspaces, 
for three different sets of initial conditions when K2 > 0.65. All orbits end at the 
attracting point P. 

Fig. 6b: The time-evolution of the Kasner solution R{t) ~ t°'^^ for the external 
space, for several values of the normalized coupling constant ^3. Again, the time 
coordinate is measured in units of 10^ tpi. Notice the slight modulation in the time- 
evolution when < K3 < 0.75. 
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